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SPECTRAL-AVERAGING ERROR REDUCTION IN 
LOW-RESOLUTION SPECTRORADIOMETRY 
by Donald R. Buchele 
Lewis Research Center 

SUMMARY 

Analytical methods are presented for reducing the error in the weighted average of 
a spectrally varying quantity when the weighting function is a spectrum measured by an 
instrmnent having low spectral resolution. The error is determined in terms of this 
resolution; this error increases as the resolution becomes lower. Spectroradiometers 
of the scanning type and the multiple -filter type are considered. 


INTRODUCTION 

This report is directed to the improvement of the weighted average of a spectrally 
varying quantity when the weighting function is a spectrum measured by use of an instru- 
ment having low spectral resolution. The instrument may have one of two forms: 

(1) A scanning spectroradiometer with a slit function of appreciable width 

(2) A filter spectroradiometer using a discrete set of fixed filters, each of appreci- 

able bandwidth 

The analysis considers analytical techniques whereby the effective resolution of each of 
these instruments may be sharpened. 

The mathematical treatment has wider applicability than just this problem. Thus, 
one may consider the weighted average Y of a spectrally varying quantity y(X) which is 
defined by 


Y = 



N(X)y(X)dX 


Jq N(X)dX 


( 1 ) 



where N(X) is a weighting function, generally called the spectrum. (Sjnnbols are de- 
fined in appendix A. ) Typically, Y may be any of the following quantities: (1) emissiv- 
ity, where y(X) is the spectral emissivity and N(X) is the spectral radiance of a black- 
body at the surface temperature; (2) absorptivity, where y(X) is the spectral absorp- 
tivity of a surface and N(X) is the spectral radiance of incident radiation; (3) trans- 
mission factor, where y(X) is a spectral transmission factor of a material and N(X) 
is the spectral radiance of the radiation entering the material; or (4) detector sensi- 
tivity, where y(X) is the spectral sensitivity and N(X) is the spectral radiance of the 
incident radiation. 

When N(X) is obtained by spectral radiometric measurements, it is known only with 
finite spectral resolution, which causes an error in the computed value of Y. It is often 
desirable to use a low-resolution spectral measurement, provided the resultant error in 
Y is small enough, because the low resolution and increased signal flux to the detector 
make possible a more favorable combination of signal-to- noise ratio and spectrum re- 
cording speed. Low resolution may be required to obtain enough signal flux for meas- 
urement. Instrument size and weight are often also reduced with low resolution. 

This report determines the relation between spectral resolution and the error in Y. 
Analytical methods are then presented for reducing the error in Y. These methods also 
improve the spectral resolution, which may be useful in other applications of spectral 
measurement. The errors in Y caused by instrumental errors such as wavelength 
calibration and flux measurement will not be considered. The error due to flux meas- 
urement is reported by reference 1 . 

The anal 3 Ttical methods for error reduction are (1) compensation of the Fourier 
transform of the spread function and (2) spectral- curve compensation. Analogous ex- 
amples of Fourier transform compensation are given in references 2 to 7. Applications 
of the technique of spectral- curve compensation are given in references 8 to 13. 

Fourier transform compensation will be treated first, because it is a general method 
applicable to both scanning spectroradiometers and filter radiometers. Then, the 
spectral- curve compensation method will be treated in its specific application to filter 
radiometers. 


GENERAL ANALYSIS 

This general analysis applies to both a scanning radiometer and a filter radiometer. 
However, a filter radiometer requires the further analysis given in the section FILTER 
RADIOMETER. 
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Spread Function 


The transmission function of a filter or a spectrometer is given by 

t(u) = t(X - X.) (2) 

where Xj is the nominal center wavelength of a filter, or a fixed wavelength setting of 
a spectrometer. Equation (2) is convenient for use in analysis; it assumes that r(v) is 
the same at all values of X^^ . Measurement of t{v) can be made in principle with a 
monochromatic source that sweeps across all wavelengths X as X. remains constant. 

As a result of finite spectral resolution, a spectral measurement N(Xj) at a wave- 
length X. is then given by a cross-correlation integral (ref. 14) of r{v) with the true 
spectrum N(X): 


N(X.) 


/ oo 

T(t')N(Xj + t>)dt' 


(3) 


(For mathematical convenience, all integrals over a wavelength will have a lower limit 
of -oo, although X cannot physically become negative. ) 

Equation (3) can be transformed into a convolution integral (ref. 14) 


N(X.) 


y* g(i^)N(Xj^ - 


(4) 


where the transmission function t(u) has been replaced by a spread function defined by 

giu) = t(-u) (5) 


The spread function is the same as the transmission function only for even functions. 
The spread function of a filter is given conveniently by equation (5) after the trans- 
mission function has been measured. The spread function of a monochromator is con- 
veniently measured by recording the spectrum of a monochromatic source fixed at X.. 

The spectrum recorded by a scanning spectrometer is given by equation (4) when 
Xj is continuously variable and replaced by X so that 


N(X) = 


g(j/)N(X - v)dv 


( 6 ) 


The spectrum measured by a filter radiometer is a set of readings, at X = X., 
given by equation (4) or (6); the readings are then interpolated to generate a spectral 
curve. Thus, the spectral curves produced by a scanning spectrometer and a filter 
radiometer agree only at X., when both instruments have the same spread function. 
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A further simplification will be applied in some of the following sections by expand- 
ing equation (6) into a sum of two integrals 

/ oo ^co 

ggMN(X - v)dv + J - p)dv (7) 

where the spread function g{p) is the sum of an even part gg(t') and an odd part g^M^ 
The odd part can be neglected, because (1) the functions generally foimd in practice are 
substantially even, and (2) the size of the odd integral can be made zero to a first-order 
approximation (ref. 2) by choosii^ the center wavelength X. of the spread fimction so 
that it is located at the centroid given by 


_ /co ^g(^)dX 

Xo 

At this value of wavelength, a moderately asymmetrical spread function jdelds very 
nearly the same measurement of N(X) as a symmetrical one. 


Averaging Error 

As a result of finite spectral resolution, the quantity computed from measurements 
is 


_ /rN(x)y(x)dx 

•XT 


Y = 


/: 


N(X)dX 


(9) 


where the result of finite -resolution measurements is designated by a bar. The error 
made in determination of Y is 


Y - 




(10) 


A useful equality (ref. 14) of the convolution integral, equations (4) and (6), (which also 
holds for the cross- correlation integral, eq. (3)) is 
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( 11 ) 


Combinii^ equations (1), 


^OO ^OO ^CO 

y.co NdX = gdi/. NdX 


(10), and (11) yields the fractional error 


AY 

Y 


/ ■OO 

^ (AN - N)y dX 


Ny dX 



where 



( 12 ) 


(13) 


Relation Between Spread Function and Averaging Error 

Substituting equation (6) into equation (12) and exchanging the order of integration 
yields 


AY 

Y 


A h(i/)g(t')di/ 
Xo N(X)y(X)dX 


- 1 


where h(t') is a cross- correlation integral 


h(i/) 


-I 


y(X)N(X - z/)dX 


The denominator appearing in equation (14) is 

h(0) = Xo N(X)y(X)dX 


(14) 


(15) 


(16) 


When g{v) is a imit-area impulse 5{v) representing infinite spectral resolution, the in- 
tegral in the numerator of equation (14) becomes h(0), and the error is zero. 

When equation (16) is substituted in (14), the error is 
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( 17 ) 


/ 'OO 

00 Mv)g(v)dv 
Y h(0) 

K g{v) is an even function and lies entirely within a strictly linear range of hCi^), the 
error is zero. 

Although both equations (12) and (17) can be used to compute the error, equation (17) 
explicitly shows the relation between the spread function (or spectral resolution) and the 
error. Equation (17) provides the means for determining how much g(v) should be nar- 
rowed, or be made more symmetrical about v = 0, should the first computation of 
I AY/y| show this quantity to be excessive. 




(b) Cross correlation of N and y. 

Figure 1. - Determination of spread-function width. 
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The spread function width that will yield an acceptably small error, yet be as wide 
as possible, can be established in the following steps (these steps will be illustrated by 
the parallel example of fig. 1): 

(1) Given a spectrum N(X) and a detector having response y(X) (e. g. , the xenon 
spectrum and the cadmium-sulphide photosurface sensitivity shown in fig. 1(a)), compute 
the cross-correlation fimction h(z^), equation (15) (fig. 1(b) for the example), 

(2) Determine a range of v about v = 0 at which the separation of h(i') from a 
tangent line at = 0 does not exceed the specified error (see fig. 1(b)). 

(3) Select a spread function half- maxi mum width that is within the selected range of 
V, and compute the error with equation (17). H the maximum error thus computed is too 
large or too small, select another width and recompute the error. 


AVERAGING-ERROR REDUCTION BY COViPENSATION OF FOURIER TRANSFORM 

Analysis 

The purposes of this section are to find the effective spread function shape produced 
by Fourier transform compensation and to illustrate a simple type of compensation. 

When an effectively even fvmction g(v) = gg(t') is assumed and the convolution symbol * is 
used, equation (6) becomes 


N(X) =gg(X)*N(X) (18) 

which has a transform 

^(s) = S>g(s)>’(s) (19) 

The shape of the spread function can be compensated to yield an effectively higher 
resolution by multiplying its transform S^g(s) with a compensating function ^q{s) so 
that ^q(s)9^(s) « constant over some range of s about s = 0. Full compensation 
(.^^(s)Sfg(s) = constant for all s) is not practical at values of s where ^g(s) « 0 be- 
cause measurement error is magnified by compensation (ref. 15). With the correspond- 
ing compensating function f^(i/), equations (18) and (19) become, respectively. 


f^(X)*N(X) =f^(X)*gg(X)*N(X) 

(20) 

jr^(s):^(s) = 

(21) 


Thus, the compensated spread function may be computed with the convolution 
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g^(X) = f^(X)*gg(X) (in order to observe the nature of the improvement). The corrected 
spectrum can be computed with the convolution f^(X)*N(X). References 4 and 5 specify 
^q(s)9^{s) for a particular application. References 6 and 7 specify .^^(s) using a 
criterion of minimum least square error of the difference N(X) - N(X), thus automati- 
cally avoiding excessive magnification of the measurement error. 

Reference 3 compensated (s) with a trigonometric polynomial which lends itself 
to simple numerical computation. It will now be applied to demonstrate the change in 
spread-function shape after compensation. For an even spread function Sq(v), the first- 
order trigonometric correction term is also an even function, and thus real, given by 

= 1 + b sin^ ffas = 1 + (b/2)(l - cos 2;ras) (22) 

where the constants b and a are constrained to make ^q(s)^^(s) « constant for small 
values of s. This constraint is satisfied by use of the moment theorem (refs. 14 and 16) 
(the second moment of the spread function is proportional to the second derivative of its 
transform at s = 0), which leads to 


ba^ = 2a^ 


where a is the standard deviation of the spread function and is given by 




^oo 

gg(z^)di^ 


(23) 


(24) 


A value of a. - o in equation (23) is suggested by reference 3, which gives b = 2. If the 
spread function has a value of a that is very large or infinite, a graphical method of 
selecting a and b must be used (ref. 3). Equation (22.) and its inverse transform to 
the V domain are shown in figure 2. 

Fourier transform compensation will now be applied, with a = a and b = 2 for the 
compensating function. Examples of Gaussian and rectangular spread functions and their 
transforms are shown in figures 3 and 4, together with the compensating function ^q(s) 
and its inverse transform. The large negative values in the compensated function ^q(s) 
of figure 4(c) may be reduced by use of a different value of a such as a = 2 (t, giving 
b = 0. 5 by equation (23). After compensation, the standard deviation of the spread fimc- 
tion always has an iindesirable increase because of the convolution However, 

the shape of the altered spread function in these examples has negative values at the 
wings, so that there is an overall improvement. When this spread function is integrated 
in equation (17), the negative values tend to reduce the error. 
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1 + (b/2)[l - cos(27ras)] 



(a) Inverse transform of compensating function. (b) Compensating function. 

Figure 2. - First-order term of trigonometric polynomial compensating function. 



(b) Inverse transform of compensa- 
ting function. 


(c) Compensating function and transforms 
of spread functions. 


Figure 3. - Example of Gaussian spread function, with standard deviation o = 0.4, compensated 
by trigonometric polynomial with constants of function ^^.(s) being a = o and b= 2. 



(a) Original and compensated spread 
functions. 



(b) Inverse transform of compensa- 
ting function. 



s 


(c) Compensating function and transforms of spread 
functions. 


Figured - Example of rectangular spread function, with standard deviation o = 0.29, compensated 
by trigonometrical polynomial with constants of function ^^.(s) being a = o and b= 2. 


Example for a Scanning Spectro radiometer 

The averaging error will be computed for the xenon-arc lamp spectrum and the cad- 
mium sulfide photosimface shown in figure 1(a). The cross correlation of these curves 
is h(t'), shown in figure 1(b). The spread function that will be assumed for the spec- 
trometer is 


gg = 1 + COS(jTI//w), I I//w| :< 1 (25) 

where w is the half- maximum width of gg(t')- This function is shown in figure 5(a). 

The error computed with equations (12) and (17) is plotted in figure 5(c). Both equations 
give the same error. 

Anal 3 Ttical reduction of the error used the compensating fvmction of equation (22). A 
value of a = a was used in equation (23) as suggested by reference 3, giving b = 2. 
Thus, the ordinates of the compensating function fQ{v) in figure 2 become -0. 5 and 2. 0 
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Spread function half-maximum width of spectroradiometer, w, jam 
(c) Averaging error (curves from eq. (12); points from eq. (N . 



Figure 5. - Scanning spectroradiometer with Fourier-transform compensation. 


with a = = 0. 36 w for the spread fxinction Convolution of tQ(v) with g^iv) 

produced the compensated fxmction gQ(v), shown in figure 5(a) for comparison with 
gg(t'). The corresponding transforms are shown in figure 5(b), and the spectra are 
shown in figure 5(d). The averaging error, plotted in figure 5(c) as a function of spectral 
resolution, has been substantially reduced. An error of 5 percent at 0. 3 micrometer 
was reduced to 0. 5 percent. Conversely, for a given 0. 5 percent error, the spread 
function could be three times wider if analytical error reduction were to be applied sub- 
sequent to the measurement. 

This example illustrates that the averaging error need not be large when the spec- 
trum consists of many emission lines rather than being a smooth ciorve. The spectral 
measurement does not have to resolve an emission line. It is generally adequate to have 
the spectral width of the spread function g(v) less than the spectral width of either NO'^) 
or y(A), whichever one has the greater width. This is shown more precisely by equa- 
tion (17), where the error becomes small when the spectral width of g(i^) is less than 
that of h(t»), and h(i^) is the cross correlation of N(X) with y(X). In the present exam- 
ple, the half- maximum width of y(A) was 0. 27 micrometer. An averaging error of 0. 5 
percent was produced with a half- maximum width g^iv) of 0. 12 micrometer, which is 
one-haK the half-width of y(X); if subsequent analytical compensation were to be used, 
the half- width of g^iv) could be 0. 38 micrometer, which is approximately equal to the 
0. 27-micrometer half-width of y(X). Thus, a small averaging error can be obtained 
with rather wide spread functions. 


FILTER RADIOMETER 

Sampling and Interpolation with a Set of Equally Spaced Filters 

The spectrum obtained with an infinite set of equally spaced symmetrical filters is 
the result of three steps: 

(1) When a spread function of a filter scans the spectrum continuously, the resultant 
measured spectrum N(X), expressed by equation (6) or (18), represents a smoothing of 
the original spectrum N(X). 

(2) K, now, the form of the spread fimction is fixed and repeated at center wave- 
lengths Xj spaced at intervals A - it represents a set of filters . Radiation 

measurement with each filter is a sample of N(X) at X - X.. The samples are 

N(X.) =N(Xq + 1A), i =-oo, . . 0, 1, . . ., CO (26) 

where Xq is the center wavelength of one of the filters. The samples are periodically 
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(d) Fourier transform of linear Interpolating func- 
tion p(v). 


Rgure 6. - Sampling and interpolation with set of filters. 


reproduced, as multiples of A are added to Xq. Figure 6(a) shows some spread func- 
tions g (u) spaced at an interval A; figure 6(b) shows the Foimier transform 9 (s) of 
one spread function. 

(3) A spectral curve is generated by fitting a curve through the sample ordinates 
N(X.). The curve-fitting method commonly used in Fourier transform analysis is the 
convolution of an interpolating function p(v) with equally spaced sample ordinates N(X.). 
The convolution integral reduces to a summation, giving a spectral curve designated by 
the subscript f: 


OO 

Nf(X) = ^ N(X.)p(i/) (27) 

i=-oo 

Figure 6(c) shows ordinates N(X.) with triangular (linear) interpolating functions p(t'). 
Figure 6(d) shows the transform ^(s) of p(^'). 

Summarizing, the three steps are (1) smoothing of N(X) by gg(i^) to produce N(X), 
(2) sampling of N(X) at intervals A to produce N(X.) and (3) interpolation of N(X.) by 
p{v) to produce Nj(X). 

The dashed vertical line in figiire 6(b) corresponds to the interval A, and its 
abscissa s^ = 1/(2 A) is the "sampling cutoff frequency. " Values of s to the right and 
to the left of s^ have different effects. 
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Values of s to the right of s. correspond to sinusoidal components, in the wave- 
length spectrum N(X), with a period less than 2A. These short-period components of 
N(X) are substantially attenuated by the inverse transform of ^(s) and then become 
components of N(X). The short-period components of N(X) are then sampled at an in- 
terval A that is greater than one-half their period; they therefore cannot be recon- 
stituted from the sampled data; they become an error of different period (called 
’’aliasing” error) in the reconstituted spectrum N£(X). The error at any wavelength 
depends on the arbitrary value of Xq in equation (26). This limitation is a result of 
the sampling theorem (ref. 14). 

Values of s to the left of s correspond to sinusoidal components, in the wave- 
length spectrum N(X), with a period greater than 2A. The long-period components of 
N(X) are attenuated by the inverse transform of 2t(s) to become components of N(X). 
These long-period components of N(X) are then sampled at an interval less than one-half 
their period; they therefore can be reconstituted in their attenuated form from sampled 
data, as proved by the sampling theorem (ref. 14). With interpolation between wave- 
lengths Xp these long-period components of N(X) can be further attenuated by the in- 
verse transform of ^(s) to become components of N^(X) (eq. (27)). Thus, the combined 
effect of the spread function and of interpolation is to attenuate the long-period compo- 
nents of N(X) by the inverse transform of product This attenuation can be 

compensated for values of j s | < s^ by either of the two general methods of reducing 
averaging error. Thereby, the long-period components of N(X) can be reconstituted. 

Attenuation by the interpolating function can be avoided by use of a sinc()Ti//A) 
interpolating function defined by 

p(t») = sinc(n-i//A) ^ sin{nv/ 1^/ {nv/ K) (28) 

which has a transform that is constant for j s j si/ (2 A) and zero elsewhere. This func- 
tion is useful in analysis, but it has practical limitations. An alternative interpolating 
function is recommended in the section Interpolation with a Finite Set of Unequally 
Spaced Filters. 

In the following sections the quantities S?g(s), <^(s), and s^ will be related to the 
spread -function shape corresponding to a set of filters and to the choice of filter spacing 
A. 


Spread Function and Fourier Transform of a Set of Filters 

The effective spread fvinction gj for each filter in a set of fixed filters differs from 
the spread function g^ for the same filter when it is used in a continuous scan mode, 
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because of sampling and interpolation. The effect of sampling is to add the ordinarily 
random aliasing error 


^(s)^g(s) for |s| > Sg 
to 

•^(s)^^g(s) for |s| s Sg 


This random aliasing error may be neglected in attempting to determine and correct for 
the averaging error. 

The effect of interpolation is to attenuate w4^(s)®g(s) for |s[ < by a factor ^(s). 
Thus, the spread fvmction for a set of fixed filters has a transform given by 


^f(s) =^g(s)^(s) ls| < Sg 

^f(s) =0 l®l > ®C 


(29) 


that includes the effect of sampling and interpolation on ^g(s). The corresponding 
spread function g£(y) for a set of filters is the convolution 

gf{v) = gg(i^)*p(i^)* sinc(7n^/A) (30) 


where the truncation at s^ is produced by sinc(7ri^/A). Equation (30) shows that g£(t<) 
is wider than g^W)- 


Width and Spacing of Filters 

The requirement that the maximum fractional error • AY/Y be less than a specified 
amovint imposes restrictions (1) on the width of the filters and (2) on the spacing of the 
filters. The recommended criteria for width and spacing are 

(1) The spread function of the filter should be the g£(i^) of equation (29) or (30), and 
its width should be selected as indicated in the section Relation Between Spread Function 
and Error. 

(2) For the filter width selected in criterion (1), the filter spacing should be close 
to twice the standard deviation a, equation (24), of the spread fimction of the filter 

gg(i^). 

The result of satisfying criterion (2) is illustrated in figvme 7, which shows rectan- 
gular, triangular, and Gaussian fvmctions with a separation A = 2 ct. These functions 
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(a) Rectangular filter function. 




(b) Triangular filter function. 



Figure 7. - Filter separation and Fourier transform of one filter. Half-maximum width of spread function = 1. 


represent filters found in practice that range from rectangular to triangular to Gaussian 
in shape. The Fourier transform ^^(s) is shown for each case. The abscissa of the 
dashed vertical line in the transform plane is the "sampling cutoff frequency, " 

Sj. = 1/(2A). This vertical line has the following significance: 

(1) Area A to the left of the line represents a correctable attenuation of long-period 
sinusoidal components in the spectrum N(X). 

(2) Area B to the right of the line represents attenuated short-period components in 
the spectrum N(X), which, after sampling, produce an uncorrectable error in the recon- 
stituted spectrum Nj(X). For the scanning radiometer, area A extends to infinity, and 
area B is zero, because the sampling interval A is zero, and s^ becomes infinite. 
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(3) The location of the line represents a compromise between several considerations: 
The correction for attenuation, area A, utilizes division by the ordinates of the curve up 
to s^; as the vertical Line is moved to the right, the distance C must not be allowed to 
become too small because measurement error for sinusoidal components of N(X) with a 
short period near to 2A would be magnified excessively (ref. 15). On the other hand, 
as the vertical line is moved to the left, area B increases rapidly, giving an undesirable 
increase in the uncorrectable error caused by sampling. 

Thus, the filter separation A must be such as to minimize the sum of variances 
represented by areas A and B. The compromise spacing A = 2 ct makes C = 0. 15 for 
the rectangular function, C = 0. 24 for the triangular function, and C = 0. 29 for the 
Gaussian function. The functions g^iv) show substantial overlapping with this spacing. 


Interpolation with a Finite Set of Unequally Spaced Filters 

The ideal band- limited interpolation with a sinc(7Ti^/A) function is not recommended 
because of the following conditions: 

(1) Filter center wavelengths are unequally spaced. Nominally equal spacing is not 
easily obtained, because of fabrication tolerances. Furthermore, the spacing is often 
progressively increased with wavelength to maintain a constant ratio of spacing to wave- 
length. 

(2) The set of filters is finite. An interpolation error arises from lack of readings 
that would be provided by nonexistent filters outside the set. These readings would con- 
tribute to interpolation within the set through the slowly decaying sine interpolating 
function. This error has been considered by references 17 and 18. Attempts have been 
made to minimize this error by the use of more rapidly decaying interpolation functions 
(refs. 19 and 20). 

(3) Extrapolation beyond the first and last filter center wavelengths is necessary to 
apply the methods of averaging- error reduction. 

(4) Negative values of radiation are not permitted. 

Cubic spline interpolation (ref. 21) is appropriate for overcoming the first three 
objections. It is described in appendix B. A function p(i>) (eq. (27)) that approximates 
cubic-spline interpolation and that is more convenient for analytical treatment of 
equally spaced intervals is 


p(i^) = 0. 4 sin(7Tt'/A)/sinh(0. 4 nv/A) (31) 

The approximation is within 1 percent everywhere except very near to Xj and to X^^. 
This function is shown in figure 8(a). For comparison, the sinc(n’t'/A) function, equa- 
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(b) Fourier transform. 

Figure 8. - Interpolating function and its transform. 


tion (28), and a linear function are also plotted. The spline fvinction has a faster decay, 
desirable for condition (2). 

The transform of equation (31), shown in figure 8(b), is 

^{s) = A sinh(2. 5 tt)/ ( cosh 2. 5 tt + cosh 5nsA) (32) 

This curve is intermediate between the linear and the sine function transforms. 

At the first and last filter center wavelengths, Xj and X^^, the cubic spline com- 
puted with equation (Bl) has zero curvature. Linear extrapolation of the slope given by 
the derivative of equation (Bl) is appropriate for condition (3). 

The negative values of radiation sometimes produced by interpolation or extrapola- 
tion are replaced by zero to satisfy condition (4). However, during averaging- error re- 
duction by the method of spectral curve compensation, to be treater later, it was found 
desirable to retain negative values for iterative adjustment of their magnitude and to 
replace negative values by zero just before integrating the product of spectral radiation 
by the spread function. 
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Example of Averaging-Error Reduction by Fourier-Transform Compensation 


Assumed conditions » - The averaging error was computed for the same xenon 
source, cadmium sulfide photosurface, cosine spread function (eq. (25)), and method of 
error reduction as used in the earlier example with a scanning spectroradiometer. The 
spread functions were spaced by the recommended 2a so that A = 2a = 0. 72 w. 

Rand om error due to sampling . - To obtain a value for the random error due to 
sampling by a given set of filters, four spectral curves were determined for that set, 
as the entire set was shifted, relative to the spectrum, by successive wavelength dis- 
placements of A/4, A/2, and 3 A/4 from an initial position. Four typical spectral 
curves are shown in figure 9 for the case of five filters having w - 0. 3 micrometer. 



Figure 9. - Xenon spectra with linear interpolation between readings of five filters of half-maximum 
width of 0. 3 micrometer, spaced by 0.216 micrometer. Each filter in set is displaced in three steps of 
0.054 micrometer to generate four curves. 


The averaging error was computed for each curve using equation (12). The average 
deviation from the mean of these four averaging errors is the random error due to 
sampling. 

Determination of compensating function . - Error reduction used the compensating 
function of equation (22), where the constants a and b were constrained by 

^^(s)^j(s) = (s) ^ constant 

for small values of s. We define 
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and choose a = cr^ as suggested by reference 3. Applying the moment theorem (ref. 14), 
we have 

= -‘®e 

4tT^O^ = -r^”(0) 

and note that 


^”(0) « -5.2 


for the spread function, 


for the cubic spline, and 


^”( 0 ) = 0 
^”( 0 ) = -2n^A^/Z 


for the linear interpolating function. Adopting the recommendation given in the section 
Width and Spacing of Filters that 

<7^ = A/2 


leads to 


for cubic spline interpolation and 


a = 0. 36 w 


a = 0. 47 w 


for linear interpolation. The second value of a provides the extra compensation re- 
quired by Linear interpolation. 

Error comparison . - The errors before and after compensation are plotted in fig- 
ures 10(a) and 11(a). An extra curve is shown in figure 10(a) for a = 0. 36 w to show 
the partial compensation when the effect of interpolation is neglected in choosing the 
value of a. Comparing these curves to the scanning spectroradiometer curves in fig- 


20 


g(v/w) Fractional averaging error, AY/Y 



Spread function half-maximum width of one filter, w, 

(a) Averaging error (curves from eq. (12) (bars are average deviation); 
points from eq. (17)). 



(b) Spread function of set of filters. (c) Transform. 


Figure 10. - Filter radiometer with linear interpolation and Fourier transform com- 
pensation. 


g(v/w) Fractional averaging error, AY/Y 


.04 — 



Spread function half-maximum width of one filter, w, pm 



Spread function half-maximum width of one fitter, w, pm 


(b) Spread function of set of filters. (c) Transform. 


Figure 11. - Filter radiometer with spline interpolation and Fourier transform com- 
pensation. 


vire 5 shows that the averaging error after compensation, for w < 0. 3, is effectively 
the same for both instruments. Figures 10 and 11 also show the random errors (average 
deviations) caused by sampling. The average deviation is doubled after correction, and 
it is then comparable to the averaging error. The average deviation is greater for 
linear interpolation than for spline interpolation. 

The spread functions and transforms, shown in figures 10(b), 10(c), 11(b), and 
11(c), are substantially different from those of figures 5(a) and (b). 


Optimum Filter Width for a Small Number of Equally Spaced Filters 

The illustrative examples have assumed that the number of filters used was so large 
that an increase in their number would not produce a significant decrease in the final 
errors. For some space applications, where mass or bulk of the instruments is of 
great significance, one may want to examine the effects of using a very small number 
of filters. This section considers this problem; specifically, the final errors associ- 
ated with using three to seven filters . 

For a given spectrvim to be analyzed, and a given nxunber of filters, there will be 
some value of w, the half- maximum filter width, which minimizes the systematic error 
after compensation, the random error due to sampling, or any chosen combination of 
these two errors, such as the sum of the systematic and random variances. This situ- 
ation is in contrast to the case where the number of filters is large, because then the 
errors always decrease as w decreases. 

For concretenesb, the problem to be treated will be confined specifically to the 
same source, detector, and filter spread function used in the previous examples and will 
assume cubic-spline interpolation. The results are plotted in figure 12 for filters that 
were symmetrically placed over the cadmium sulfide response band of 0. 45 to 0. 9 mi- 
crometer. Using the criterion of minimum sum of systematic and random variances, 
one result is that the total variance is a minimum, or close to a minimum, for 
w » 0.25 micrometer, regardless of the number of filters. Another result is that, on 
the locus of minimum variance, four filters have a wavelength span 50 percent greater 
than the 0. 45-micrometer span of the detector. As the number of filters increases, 
the wavelength span of the filters approaches that of the detector. In this example, 
there is little advantage to using more than four filters. 
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AVERAGING-ERROR REDUCTION BY SPECTRAL-CURVE COMPENSATION 


Analysis 

After readings have been taken with a set of filters, a low- resolution spectral curve 
is drawn by interpolation between the readings when these have been plotted at the center 
wavelengths of the filters. Using this low-resolution spectral curve in place of the real 
spectrum, readings of the filters may be calculated. It is found that the measured and 
the calculated readings do not agree, as they would if the low-resolution spectrum were 
equivalent to the real high-resolution spectrum in the production of filter readings. 
However, references 8 to 13 obtain such equivalence by applying corrections to the 
measured readings in such manner that the corrected spectral curve yields calculated 
readings equal to the measured readings. This corrected low-resolution spectrum is 
named the ’’filter equivalent spectrum” by reference 8. An ’’effective spread function” 
gg(t') associated with the equivalent spectrum method will be found, and it will be shown 
that the transforms of both the filter spread function and the interpolating function are 
compensated in the range |s| < s^. Thereby, the averaging error is reduced. 

The measured filter readings are a sequence of impulses that have a transform 
J/ (s) (ref. 16) which is periodic, with a period A, and is given by 

OO 

^*(s) = ^ J"[s + (n/A)]^g[s + (n/A)J (33) 

n=-oo 


The equivalent spectrum is generated by fitting a curve through the corrected ordi- 
nates spaced at intervals A. Curve-fitting with an interpolating function p(t') and the 
corrected ordinates Ng(X^) generates the equivalent spectrum 


Ng(X) - ^ ^ Ng(Xj^)p(i^) (34) 

i^-co 

The calculated filter readings are the result of integrating Nj,(X) • gg(t') to produce 
readings Ng(xp. The samples are a sequence of impulses that have a transform like 
equation (33), given by 


CXD 

./f'g(s) ./Kg[s + (n/A)]^g[s + (n/A)] 

n=-oo 


(35) 
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The method of spectral curve compensation makes the filter readings equal to each 
other, and thus their transforms become equal, so that 

^*(s)=>*(s) (36) 

K equations (33) and (35) are equated to each other, ^^(s) cannot be canceled if the terms 
of the summation overlap, which they generally do. In general, ^(s). 

The effective spread function for the equivalent spectrum method has a transform 
defined by 

^g(s) s ./<^g(s)/J^(s), |sl<s^, (37) 

In general, ^g(s) 1, because equations (33) to (36) are complex and asymmetrical. 

Furthermore, ^g(s) depends on the spectrum as well as on the filter spread function. 
However, an approximation of ^g(s) can be made, if the spectra N(X) and Ng(X) are 
sufficiently smooth curves, so that the terms in equations (33) and (35) do not signifi- 
cantly overlap. Then equation (37) becomes ^j-(s) = 1. Then the spread function 
gg(t») = sinc(7Tt'/A). Thus, the filter spread function gg(i^) has been fully compensated 
in the range |s| < s^. Since interpolation preceded the generation of Ng(X) and 
(eqs. (34) and (36)) and compensation followed the generation of these quantities, no 
separate compensation of the interpolating function is required. 

The amount of compensation applied by this method must be limited by the proper 
choice of filter spacing (A = 2a ) so that the length C in figure 7 is not too small, as 
was discussed in the section Width and Spacing of Filters. 


Example for a Filter Radiometer 

The example used is the same as in the section Example of Averaging Error Reduc- 
tion by Fourier- Transform Compensation. The Newton-Raphson method was used to 
adjust the original filter measurements to provide the filter equivalent spectrum and 
thereby to yield zero differences between the originally measured filter readings and the 
filter readings calculated for the filter equivalent spectrum. Fast convergence of the 
method was obtained by use of a tridiagonal Jacobian matrix, with the elements computed 
by numerical differentiation. Convergence was rapid because the size of the partial 
derivatives falls rapidly away from the principal diagonal. The subroutine of appendix B 
was used to solve the Jacobian system of equations. 

The error is plotted in figure 13 for both linear and spline interpolation. The un- 
compensated values shown are the same as those shown in figures 10(a) and 11(a). With 


25 


. 08 , 




Spread function half-maximum width of one filter, w, pm 

(a) Linear interpolation. ib) Spline interpolation. 

Figure 13. - Filter radiometer averaging error with spectral curve compensation. (Curves from eq. (12) (bars are average 
deviation); points from eq. (17), compensation with g^(v) = sinc(7rv/A), eq. (28).) 


compensation, the points computed by equation (17) show overcompensation relative to 
the curve computed by equation (12) for w > 0. 2 micrometer. Thus, the correct spread 
function for the equivalent spectrum probably lies between the function sinc(7rt'/A) and 
the spread functions in figures 10(b) and 11(b) for Fourier transform compensation. 

Both systematic averaging error and random average deviation, after spectral- 
curve compensation, are less than the corresponding errors after Fourier transform 
compensation that are shown in figures 10(a) and 11(a). At 0.25 micrometer, a 
3-percent systematic error was reduced to 0.3 percent, and the average deviation was 
0.3 percent. Both linear and spline interpolation jdeld the same averaging error, but 
the average deviation is smaller with spline interpolation. Thus, the smallest errors 
are obtained with spline interpolation and spectral-cmve compensation. 

Figure 14 shows how the ordinates of the filter equivalent spectrum differ from 
those of the faired spectral curve drawn through the original measurements. The com- 
parison is shown for both linear interpolation (fig. 14(a)) and cubic-spline interpolation 
(fig. 14(b)). A large negative value of the ordinate for the first filter, at X = 0. 175 mi- 
crometer, was required to generate the filter equivalent spectrum. As indicated in the 
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Compensated 




(b) Cubic spline interpolation. 

gure 14. - Comparison of linear and spline interpolation with five filters of half-maximum width of 0.3 micro- 
meter and xenon spectrum. 



section Interpolation With a Finite Set of Unequally Spaced Filters, the negative values 
were replaced by zero before integrating the product of spectral radiation and the spread 
function. 


SUMMARY OF RESULTS 

An equation was derived for the error of a spectral average caused by low spectral 
resolution. In this equation, the spread function appeared explicitly so that the error 
could be related to the spread-function width. Error reduction by Fourier-transform 
compensation and by spectral- curve compensation was analyzed and examples were 
given. Both methods gave similar results. In an example of spectral- curve compensa- 
tion, a systematic averaging error of 3 percent was reduced to 0. 3 percent, with a 
random average deviation of 0. 3 percent. Spectral- curve compensation is more adapt- 
able than Fourier-transform compensation in those practical applications where the 
spread functions are not alike or are not equally spaced. Error reduction is also a use- 
ful approach to estimating the averaging error. The procedure is to find the value of 
the average before and after appl 3 n.ng error reduction. Trial examples show that the 
error is less than the difference between these averages. 

Spectroradiometers of the scanning type and the multiple-filter type were con- 
sidered. Criteria were established for the optimum spacing of filters, having a given 
spread- function width, to minimize the errors. Recommendations were given for cal- 
culating the optimum location of the center wavelength of a filter and for interpolating 
between filter readings. An example showed that (1) spline interpolation produced a 
smaller averaging error than linear interpolation before averaging- error reduction was 
applied, (2) linear and spline interpolation produced substantially equal averaging errors 
after averaging- err or reduction was applied, and (3) spline interpolation produced a 
smaller average deviation of the averaging error both before and after averaging- error 
reduction was applied. 

The example also showed that, when a very small number of filters is used, there 
exists an optimum spread-function width that minimizes the averaging error or the over- 
all error. In the particular example chosen, there was also little advantage to the use of 
more than four filters. 

Thus, using the optimum filter- selection criteria and techniques of data analysis, 
small overall errors are achievable even when a small number of low-resolution filters 
is used. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 29, 1970, 

720-03. 
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APPENDIX A 


SYMBOLS 


A normalizing constant, eq, (13) 

a constant of function ^^(s), eq. (22) 

b constant of function eq. (22) 

C amplitude of Fourier transform at s = 1/ (2A) 
d filter spacing for spline interpolation, eq. (B2) 

^ Fourier transform of f 

f fimction to cojnpensate g 

^ Fourier transform of g 

g spread function 

h cross correlation of N and g, eq. (15) 

N infinite resolution spectrum 

^ Fourier transform of N 

^ Fourier transform of p 

p interpolating fimction 

s abscissa of Fourier transform 

u variable for spline interpolation 

w half- maximum width of spread function 

Y weighted average of y 

AY error in Y 

y response function 

z variable for spline interpolation 

A separation of filter center wavelengths 

X wavelength 

X. nominal center wavelength of a filter, or wavelength setting of a spectrometer 

V wavelength auxiliary variable, = X - X. 

a standard deviation of spread function, eq. (24) 

T transmission factor, or function t(v) 


29 



Subscripts: 

C compensated 
c samplii^ cutoff 
E equivalent spectrum 
e even spread fxmction 
f set of filters 

0 odd function 

p interpolating function 
0 reference value 

Superscripts: 

~ associated with low resolution 
* Fourier transform of a sequence of impulses 
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APPENDIX B 


CUBIC-SPLINE INTERPOLATION 


The cubic spline is a cubic polynomial between each pair of ordinates, which also 
has a continuous second derivative that varies linearly between each pair of ordinates. 

It is computed within each interval between filter readings N. and N.^j at center 
wavelengths Xj and X.^j^ by the equation 

N(X) = [(> - + (Xj^, - 

- (X - Xj)(Xj^j - x)[(dj + Xj^j - X)uj + (d. + X - Xj)uj^j/edj I ' 

' Li = Ij • • • , n - 1 

(Bl) 


The computation is made for n filters with spacings 




i = 1, . . . , n - 1 J 


(B2) 


The second derivatives Uj of the spline are the solution of the system of equations 


^2 ^2 ^ ^ 


^^2 


^2 

■^3 ®3 ^3 ® 


"3 


^3 

1 

o o 

o ! 

: 

: 

1 

1 


%-l_ 


• • • 
Vi, 


(B3) 


with = Uj^ = 0 and 


Zi = - N,)/di - 6(5jj - N,.i)/dj., 1 = 2, . . . , n - 1 

The three diagonals of the matrix are 


(B4) 


Bi = 2(dj_i + dj) 

C. = d. 


> i = 2, . . . , n - 1 
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An efficient method for the solution of a tridiagonal system Uke equation (B3) is 
given by reference 22. The solution is concisely expressed as the following FORTRAN 
IV subroutine: 


SUBROUTINE SIMEQ (NtA.UJ 

C SOLVES TRIOIAGONAL N BY N MATRIX, MAIN DIAGONAL A( 2, I » , 1= 1, N, LOWER DIAGONAL 
C A(1,I),I = 2,M , UPPER DIAGONAL A < 3 , I ) ,I = 1 ,N-1 , RIGHT SIDE A( A , I) , 1= 1, N ) , 

C SOLUTION U(I),I = 1,N, STORAGE A (5 , I) ,A ( 6 , 1 ) , I = I , N 
DIMENSION A(6,N), U(N» 

A( 6, 1 )= A ( 2, 1) 

DO 10 1=2, N 

A( 5, n= A ( I, I )/A(6,I-l » 

1C A( 6, I ) =A( 2,1 )-A (5,1 )*A (3,1-1) 

U( 1 )=A(4, 1) 

DO 11 1=2 ,N 

11 U(I)=A(4,n- A(5,I)*U(I-1) 

U(N )= U(N )/A(6,N) 

00 13 J=2,N 
1= N-J+1 

13 U( 1 ) = ( U( I )-A( 3,1 )*U(I+1) > /A(6,n 
RETURN 
END 


To solve equation (B3) for u^^ using equation (B4), let N = n - 2 and 

ua)^u.^^ 

A(1,I) =d. 

A(2,I) = 2(d. + d.^j) 

A(3,I) ^d.^1 
A (4, 1) = 

where I = i = 1, . . . , n - 2. 
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